Bond Percolation in Higher Dimensions 
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We collect together results for bond percolation on various lattices from two to fourteen dimensions 
which, in the limit of large dimension d or number of neighbors z, smoothly approach a randomly 
diluted Erdos-Renyi graph. We include new results on bond diluted hyper-sphere packs in up to 
nine dimensions, which show the mean coordination, excess kurtosis and skewness evolving smoothly 
with dimension towards the Erdos-Renyi limit. 
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I. INTRODUCTION 

Percolation theory [l|, Q asks if there is a connected 
path across a system. Examples are water percolating 
through ground coffee beans and forest fires spreading 
from tree to tree. In order to have a control parame- 
ter, percolation is often studied on lattices of dimension 
d greater than 1, where percolation disappears when the 
random removal of bonds has decreased the bond concen- 
tration p to a critical value Pc- The disappearance of per- 
colation is well studied and is a second order phase tran- 
sition with a set of critical exponents that obey scaling 
laws ^^. Our interest here is in the values of Pc studied 
over diverse geometries. Although this is an old subject, 
interest continues, including in higher dimensions, where 
rigorous bounds on pc have recently been established [3] . 
We denote by z the number of initial bonds at any site 
of a particular regular lattice, (e.g. triangular net, sim- 
ple cubic, etc.) before bond dilution occurs {p = 1). It 
is convenient to define the mean coordination (r) at the 
percolation point as 



ZPc 



(1) 



which facilitates the comparison between various lattices 
in various dimensions as the mean coordination (r) at 
percolation varies much less than pc itself. A very simple 
argument suggests that (r) = 2 at the transition as each 
site must have one bond entering and one bond leaving 
to form a connected pathway. While this is the most 
efficient scenario, it does not happen quite this way in 
a random system for two reasons. First, there is redun- 
dancy where there is more than one connection between 
two points, leading to a loop. Loops push the mean co- 
ordination (r) above 2, because at least two sites with 
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coordination 3 must be involved in forming a loop. Sec- 
ond, there is irrelevancy where dead-ends and isolated 
regions are formed that would not carry a current if the 
bonds were wires in a conducting network. Irrelevancy 
pushes the mean coordination (r) below 2, as some sites 
are singly coordinated. Both these situations are illus- 
trated in Figure [TJ We will see that there is a tendency 
for the effects of redundancy and irrelevancy to cancel 
making (r) = 2 a not unreasonable starting approxima- 
tion for low dimensions d and/or low initial coordination 
z. However, for very high dimension d or initial coordi- 
nation z, the mean coordination number (r) approaches 
unity because of the preponderance of dangling bonds. 

The result (r) = 2 at the transition can also be derived 
by Maxwell type constraint counting [4| of the number 
of floppy modes [5|, Q or residual degrees of freedom / 
in the system. Connectivity percolation, which is the 
subject of this note, can be regarded as a special case 
of a larger class of problems where instead of having a 
single degree of freedom per site there are g degrees of 
freedom. An example would be vector displacements in 
two dimensions, where g = 2 . For g > 2 suchproblems 
are usually referred to as rigidity percolation [5|, |6[ . 

Maxwell constraint counting [J| is more usually em- 
ployed in problems involving rigidity, but can also be 
applied to connectivity percolation problems as a special 
case, with g = 1. More generally, there are g degrees of 
freedom associated with each site and z constraints are 
present (the number of bonds at each site is assumed to 
be exactly z everywhere initially) with probability p, so 
that 



f = g- zp/2, 



(2) 
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which goes to zero at pc = 2g/z, and hence gives the 
result (r) = 2g at percolation. Note that the number 
of floppy modes is not exactly zero at the transition as 
fluctuations in local coordination number allow for local 
redundancy and irrelevancy, but nevertheless it has been 
shown that the number of floppy modes at the transi- 
tion is extremely low [7|, making (r) = 2g an unusually 



accurate approximation for g > 2 (typically within one 
percent). For example, in the case of rigidity percolation 
of a triangular net under bond dilution Maxwell count- 
ing gives a result of (r) = 4 while numerical simulations 
find (r) = 3.961 ± 0.002, which is very close to, but 
clearly less than, 4. However, the constraint counting 
result (r) = 2 for connectivity percolation gets worse in 
higher dimensions in which (r) = 1 is reached. Neverthe- 
less, (r) is a more useful variable than pc as it changes 
much less rapidly with dimension, and we will focus on 
it here. 




FIG. 1. Showing a connected path in black across part of 
a sample with redundancy via the red loop which is over- 
constrained with one redundant bond and irrelevancy via the 
blue region that contains dangling ends that are not involved 
in percolation. 



II. BETHE LATTICE 

A useful universal guideline is provided by the Bethe 
lattice which is a tree-like network that contains no loops 
as illustrated in Figured] @: Q • If each node of the tree is 
z coordinated before dilution, then for a connected path 
there must be one way in from a previous layer, and one 
of the remaining 2 — 1 ways out must be occupied, so that 
Pc = l/{z — l); a result which can be rigorously found |9[. 
Hence the mean coordination (r) at percolation is given 
by 



(r) 



ZPc 



z/(z-l). 



(3) 




FIG. 2. Showing a tree or Bethe lattice, reproduced from 
reference Hi 



Of course not all sites have exactly this coordination 
as there is a binomial distribution of local coordination 
numbers due to the random dilution; so the probability 
of a site having r bonds present out of a total of z possible 
is given by 



Pir) 



E 

r=0 



p^ii-py- 



(4) 



and hence the n moment (r^) is given by 

z 

leading to the mean coordination 

(r) = zp 
and the square of the width Ar given by 

{/\rf = {r^)-{rf = zp{l-p) 



(5) 

(6) 

(7) 



This expression for the width is quite general for any 
network with fixed initial coordination z at every site, 
upon random bond dilution. For the Bethe lattices at 
the percolation threshold, this width becomes 



Ar 



^z{z-2) 
(z-1) 



(8) 



A particularly interesting limit is large z ^ oo where we 
obtain what we will refer to as the Erdos-Renyi limit; 
reached when percolation occurs upon bond dilution in a 
graph that initially has N nodes, each one connected to 
every other node [? ] as iV — >■ oo. In this limit 



= Ar = 1. 



(9) 



This is the limit of a large graph of nodes, where every 
node is connected to every other node with probability 
p. In the limit that the number of nodes goes to infinity, 
the chance of finding a loop becomes infinitesimally small 
and hence the large z Bethe lattice result is obtained. An 
example of a finite Erdds-Rcnyi graph [? ] is shown in 
Figure @. 



III. A UNIVERSAL PLOT 

Using what we have jotted down in the previous para- 
graph, it is convenient to combine all results for bond 
percolation on various lattices as a plot of the mean co- 
ordination (r) against the width of the distribution, or 
variance, Ar which is shown in Figure ((31). 

The results for the 2d, 3d and hypercubic lattices 
are conveniently summarized with original references in 
Wikipedia |10j . The two dimensional results, shown in 
red in Figure (|3]), are from left to right, following the thin 
red line, honeycomb, kagomc, square net and triangular 




FIG. 3. Showing a bond diluted finite Erdos-Renyi graph, 
where before dilution every node was connected to every other 
node. 



net. The general trend is higher initial coordination z to 
the right going to lower initial coordination z to the left, 
which tends to the isostatic point shown at (0,2). The 
point for the kagome lies above the point for the square 
net in the center and gives an idea of the (modest) effect 
of the detailed lattice structure as both have sites with 
four neighbors initially. Nevertheless the overall trend 
that the red points get closer to the isostatic point as the 
initial coordination z is decreased is clear. 

The three dimensional results, shown in blue, are from 
left to right, following the thin blue line, diamond, sim- 
ple cubic, body centered cubic and face centered cubic; 
with the latter two close together but following the gen- 
eral trend with higher initial coordination z to the right 
going to lower initial coordination z to the left, which 
again tends to the isostatic point at (0,2). Also included 
in Figure (JH) are the results for bond diluted hypercubic 
lattices from d = 2 up to d = 13 where the mean coor- 
dination (r) is obtained from ([TJ and the variance from 
0. 

The results for diluted non-crystalline hypcrsphere 
packings are new and were obtained from computer 
simulations of jammed configurations of iV = 262144 
monodisperse particles (in 2d a 50-50 mixture of bidis- 
perse particles with size ratio 1.4:1 was used to avoid 
crystallization) as described in reference [11|. The par- 
ticles interact with a harmonic contact potential defined 
as 



V{r) =e(cr-r)^e(fT-r) 



(10) 



where a is the particle radius, e the energy scale of the 
potential, and r the distance between particles. En- 
ergy is minimized at a given packing fraction via either 
a conjugate grad ient [12| or fast incrtial relaxation en- 
gine (FIRE) |13| minimization technique. Starting from 
a random configuration at a density well above jamming 
and given two values of packing fraction that bracket the 
jamming transition density the jamming point is found 
via a golden mean bisection search. Jamming is identi- 
fied as the packing fraction corresponding to the onset of 
non-zero energy as derived from the potential (jlO[) . For 



the purposes of percolation studies, two hypersphcres are 
said to be connected neighbors if there is a non-zero over- 
lap between them. We find that the values of (r) and Ar 
are rather insensitive to a (small) distance from the jam- 
ming transition. 
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FIG. 4. Showing results for the mean coordination against the 
variance for 2d lattices (red), 3d lattices (blue), hypercubic 
lattices (black) and random hypersphere packings (gray), at 
the percolation threshold. The straight lines joining adjacent 
points are only for guidance of the eye. The green line is the 
Bethe lattice result with the isostatic point at (0,2) and the 
Erdos-Renyi result (1,1) shown as the large purple dot. The 
dashed line shows the result of a 1/(2 — 1) expansion [Ij given 
in equation pip for hypercubic lattices 

Note that both sets of high-dimensional results, for 
bond percolation in hypercubic lattices and random hy- 
persphere packings approach the Erdos-Renyi limit, as 
can be seen from Figure (|3]). The highest dimension ex- 
plored of d = 13 for hypercubic lattices and d = 9 for 
random hypersphere packings are already very close to 
the point (1,1). This is because loops become less impor- 
tant in higher dimensions, discussed next. 

If we consider the mean and variance of the percolation 
variable r then the Erdos-Renyi limit is (from Eqns. [3]-|ni) 
the same as the z-going-to-infinity limit of the Bethe lat- 
tice (tree) result. The tree is, in turn, the looplcss limit of 
a general lattice; and, from simple geometric path count- 
ing considerations, the looplcss limit is the large z, and 
equivalently the large d, limit of a general lattice. The 
probability of two sites being joined by a graph with n 
links will be proportional to p". Now consider all graphs 
with n steps. For the "trees" we have r = n/2, and for 
all other graphs with a partial or full loop r > n/2. The 
key observation is that in high dimensions, p goes like 



1/d, as can be seen for the Bethe lattice in Equation (8). 
For example, hypercubic lattices have z = 2d and for 
random packings z is even larger. Therefore as d goes 
to infinity those diagrams with r = n/2 overwhelmingly 
dominate and hence only the trees contribute, and the 
Erdos-Renyi limit is reached. While this is not a formal 
proof, it demonstrates the plausibility of the result, and 
should form the basis for a formal mathematical proof. 

For completeness, we include the results of Gaunt and 
Ruskin [IJI who performed a l/(z— 1) expansion for bond 
percolation on bond diluted hypercubic lattices where 
z = 2d and found that percolation occurs at 



Pc = a[l + -a^ 



15 



57a'' 



(11) 



where cr=(2; — 1)~^. Note that the leading term is the 
Bethe lattice result. We include this result in Figure (j4]) 
as a dashed line, which is seen to be very close indeed 
to the results of numerical simulations (black dots) for 
hypercubic lattices with d > 3, then deviating at d = 2 
for the square lattice. 

Another convenient way to monitor the approach of 
dilute hypercubic lattices to the Erdos-Rcnyi limit, is to 
track the skewness 71 and excess kurtosis 72 which re- 
spectively monitor the evolution of the asymmetry and 
the deviation from Gaussian behavior of the distribution 
of contacts (for a Gaussian distribution 71 = 72 = 0). 
These are defined in terms of the moments of the distri- 
bution as 



7i 



72 



{{r-{r)f) 
((r- (r)) ) 



((r-(r))^) 



((r-(r))^) 

For Bethe lattices they take the values 

(3-2(r)) 
VW(2-(r)) 



7i 



1 



72 



(r)(2-(r)) 



6 ((0-1) 
(r) 



(12) 



(13) 



(14) 



(15) 



These are plotted as the solid lines in Figure ([5]). In the 
limit of a Bethe lattice with large z, the distribution of co- 
ordination number becomes a Poisson distribution with 
p(r) = e~^/r! and thus (r) ~ Ar = 71 = 72 = 1- Note 
that for the Bethe lattice, the skewness goes through zero 
at (r) = 3/2 which corresponds to z =3, and the excess 
kurtosis goes through zeros at (r) = (9 ± \PS) 6 = 1.211 
and 1.789 which corresponds to z = 4 ± \/3 = 2.227 and 
6.928 respectively. 

The skewness and the excess kurtosis for the hypercu- 
bic lattices can be obtained for the known values of Pc 
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FIG. 5. Showing the skewness (red line) and the excess kur- 
tosis (green line) as a function of the mean coordination (r) 
for Bethe lattices at the percolation threshold. Also shown 
are the skewness (triangles) and excess kurtosis (squares) for 
hypercubic lattices as gray symbols and random hypersphere 
packings as black symbols. The straight lines joining adjacent 
symbols are guides to the eye. The Erdos-Renyi result (1,1) 
is shown as the large purple dot. The dashed line shows the 
result of a 1/(2 — 1) expansion [iJl given in equation (|lll) 



from reference [1^ and using equations (|4]) and (O re- 
spectively. For the binomial distribution, the skewness 
is 



71 



l-2p 



-Jzv (1 - v) 



and the excess kurtosis is 



72 



6 



zp (1 — p) z 



(16) 



(17) 



and these are also plotted at the percolation threshold in 
Figure ([S]), which shows how they approach the Erdos- 
Renyi limit in high dimensions, providing further evi- 
dence of the relative unimportance of loops in connec- 
tivity percolation in higher dimensions. Results for the 
skewness and excess kurtosis can also be obtained from 
the expansion [ij] given in equation pTj) . coupled with 
equations p^ and p^ . and are shown as the dashed 
lines in Figure ([S]). Also shown in Figure [5] arc directly 
computed results for the skewness and excess kurtosis 



d 


(r> 


Ar 


Skewness 


Excess Kurtosis 


2 


1.9174 


1.3373 


0.1687 


-0.3890 


3 


1.4435 


1.2274 


0.6113 


0.1368 


4 


1.2289 


1.1234 


0.7749 


0.4499 


5 


1.1338 


1.0682 


0.8535 


0.6242 


6 


1.0890 


1.0423 


0.8954 


0.7244 


7 


1.0642 


1.0292 


0.9214 


0.7891 


8 


1.0459 


1.0181 


0.9397 


0.8355 


9 


1.0294 


1.0104 


0.9596 


0.8877 



TABLE I. Tabulated values for (r), Ar, Skewness, and Excess 
Kurtosis for random hypersphere packings of A'^ = 262144 
particles in dimensions d = 2 — 9. Note that all packings are 
constructed with monodisperse spheres except for d = 2 for 
which a 50-50 mixture of bidisperse particles with size ratio 
1.4:1 is used. 



for bond-diluted random hypersphere packs at the per- 
colation threshold in higher dimensions. Again a similar 
trend towards the Erdos-Renyi limit in high dimensions is 
very apparent. All results for bond-diluted hypersphere 
packings at the percolation threshold are tabulated in 
Table dH. 

These kinds of argument extend from percolation to a 
range of other processes. Among them are other q state 
Potts models (the g — >■ 1 limit is bond percolation [15|). 
which includes the Ising model [q = 2). This was per- 
haps the first system for which small \Jz was systemati- 
cally exploited by Brout and Englert [la, [l3l • The limit 
l/z — ?► gives mean field theory, associated with the tree 
graphs of the linked cluster many-body theory. This is 
the starting point for a 1/ z expansion involving graphs 
with increasing numbers of loops, which account for the 
fluctuation effects absent from mean field theory. It is 



interesting to note that similar arguments to those given 
here were previously given by Brout [lg| who exploited 
the link between trees and mean field theory for the Ising 
model, using large z, where the factor {JksT)" in an n*'' 
order graph being analogous to the pc" here. In the Ising 
model J is the exchange interaction between spins and T 
is the temperature. 

The related role of higher dimensions reducing fluctu- 
ations is of course well known in such contexts [3, [ij] , 
as is its role in reducing the probability of returns to the 
origin (loops) in random walks and related dynamic pro- 
cesses. These aspects suggest future work exploiting the 
approach used here for other processes. 



IV. CONCLUSIONS 

We have shown that all bond dilution results have uni- 
versal features so that results for various lattices in var- 
ious dimensions can be displayed on a single plot and 
these results approach the Erdos-Renyi limit in high di- 
mensions. The Erdos-Renyi limit is when percolation 
occurs upon bond dilution in a graph that initially has 
A^ nodes each one connected to every other one [? ] as 
N ^S' oo. It is also shown here that the mean coordina- 
tion at percolation (r) is often a more useful universal 
parameter than the percolation concentration itself pc- 
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